One- Dimensional Disordered Supersymmetric 
Quantum Mechanics: A Brief Survey 
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1 Introduction 

Many physical systems can be described using the concept of a random Hamil- 
tonian. Such a formulation is useful when the Hamiltonian depends on a set of 
quenched variables. In most cases it is convenient to describe them by a set of 
random variables distributed according to some probability law. Consider, for 
instance, the quantum dynamics of a particle interacting with randomly dis- 
tributed scatterers. If the potential is a sum of two-body potentials the Hamil- 
tonian may be taken in the form 



Here the quenched variables are the positions of the impurities and the number 
N of scatterers. They can, for instance, be modelled by independently random 
variables distributed according to a Poisson distribution. 

In this context one of the most elementary quantities of physical interest is 
the average density of states p{E)- If the potential is repulsive and short range 
then p{E) vanishes exponentially at the bottom of the spectrum 




p{E) 




E-»0 
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where d is the dimension of the space. This non-analytic behaviour was first 
discussed by Lifshits (Lifshits 1965) and then studied by a number of authors 
(see for example Luttinger and Waxier 1987). The physical mechanism which 
leads to this behaviour is the occurence of large regions of space that arc free 
from impurities and where the particle can move freely. Although these are 
exponentially rare events, they nevertheless contribute in the thermodynamic 
limit. This singular behaviour may also be derived using instanton techniques 
(Neuberger 1982). 

However there do exist some systems for which the density of states has a 
very different behaviour: 

1. The vibrations of a chain consisting of harmonic strings and random masses 
gives a spectral density with an accumulation of states at low frequency. 
This model, first introduced by Dyson (Dyson 1953), is in fact equivalent to 
the onc-dimcnsional Anderson model with off-diagonal disorder. 

2. In particle physics, the investigation of the random Dirac operator has stim- 
ulated interesting conjectures related to chiral symmetry breaking (Floratos 
et al. 1980). A particular model in dimension 2+1 is to take as a Hamiltonian 
the square of the Euclidean Dirac operator coupled to a random magnetic 
field H = = ~{df^+iAf^)'^ + -^a^i,F'^'' . In this case the low energy density 
of states must fullfil the inequality (Casher ct al. 1984) p{E) > Po{E), where 
Po{E) is the free density of states. This means that there is an accumulation 
of low energy states which contributes to the chiral condensate. 

3. The ono-dimensional Schrodinger Hamiltonian H = + <f>'^{x) + a^cj)' {x) 
which was first introduced by Witten (Witten 1981) as a toy model of su- 
persymmetric quantum mechanics (for a review see for example Cooper et 
al. 1995) provides a localization model with very unusual low energy spcic- 
tral properties. In certain cases it gives rise to an accumulation of levels 
at low energy. The density of states displays either a power law behaviour 
p{E) ~ E'^~^ or a logarithmic singularity p{E) ~ ^^^y^ of the same 

form as in the Dyson model (Dyson 1953). 

A common feature of these three models is the fact that the zero energy 
wave function is exactly known for any realization of the disorder. Using this 
property we have presented, for the supersymmetric model (Comtet et al. 1995), 
a physical picture that accounts for the different behaviours of the density of 
states at the bottom of the spectrum. We believe that this model is probably 
generic, by which we mean representative of a whole class of systems in which 
the disorder is encoded in the ground state. Since this model is easier to handle, 
because one can use a wealth of techniques specific to one-dimensional systems, 
we will concentrate on this case. The recent literature shows a revival of interest 
for these problems, mainly in the context of condensed matter physics. We will 
briefly comment on this work and also draw attention to earlier work which is 
scattered in the literature and has so far remained unnoticed. 
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In part 2 we will review the basic mechanism which leads to these singularities 
and underline the differences with the usual Lifshits singularities. We will also 
mention some recent applications of the supersymmetric model to quantum spin 
chains. Applications to classical diffusion in a random medium (Bouchaud et al. 
1990, Bouchaud and Georges 1990, Oshanin et al. 1993b) will not be discussed 
here. 

In part 3, following (Brodcrix and Krcc 1995) and (Shclton and Tsvelik 1997) 
we compute the correlation function of the zero energy states. 

In part 4 we discuss the fluctuation properties of the ground state wave func- 
tion using two different approaches. When ^(x) is white noise the wave func- 
tion is an exponential functional of the Brownian motion. Such functionals were 
studied extensively both in the mathematical (Yor 1992) and physical literature 
(Monthus and Comtct 1994, Comtct ct al. 1996). We will use our previous work 
to compute their statistical properties. 



2 Spectral properties 

The one-dimensional Schrodinger Hamiltonians 

H± = -^ + ^\x)±<l>'{x) (1) 

may be rewritten in the factorizcd forms = Q^Q and = QQ\ where 
Q = —-^ + (j>{x). This implies that and H- have the same spectrum for 
E > Q. When 4>{x) is random, they are characterized by the same localization 
length and density of states. These quantities have been computed exactly in 
two cases for which we now recall the main results. 



2.1 White noise potential 

{'P{x)) = Fo 



{(t){x)(t){x')) - {(l){x)f =aS{x- x') 

The integrated density of states N{E) and the localization length X{E) are 
respectively 

A-'(£) = -f + 

where z = ^ and ^ = J^(z) and N^{z) are Bessel functions. 

Equation (2) was first obtained by Ovchinnikov and Erikmann (Ovchinnikov 
and Erikmann 1977) and then rediscovered independently in (Bouchaud et al. 
1987). These results can be derived either by the node counting method (Luck 
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1992) or by the replica trick (Bouchaud et al. 1990). By this latter approach one 
can also compute the Green's function at non- coinciding points. 

The low energy behaviour of the density of states and the localization length 
are given for = by 

N{E) ^ 

X{E) ~ -InE 

and by 

N{E) ~ Ei^ 
X{E) ~ E'^-i 

E^O 

when /U ^ 0. 



2.2 Random telegraph process 

The function (^(x) is described by an ensemble of rectangular barriers (Bender- 
skii and Pastur 1974) with alternating heights ipo and ipi of random length / 
distributed according to an exponential law Po,i(0 — no,ie~"°'i'. 

In the case (f)o = —<f)\ this model yields the same low energy behaviour as 
above (Comtet et al. 1995). The parameter /i is now given by /i = . The 

main interest of this model is to provide a physical picture of the low energy 
behaviour. 

For /i 7^ the potential V{x) = (fy^ + (f)' is constant everywhere except at 
the positions where (j){x) has a discontinuity. One thus obtains a sequence of S 
functions with alternating signs. The attractive 5 potentials can support bound 
states which would have exactly zero energy if one would ignore the couplings 
to the other peaks. By taking carefully into account these couplings, one can 
recover the low energy power law behaviour (Comtet et al. 1995). The physical 
picture that emerges from this analysis is that the low energy states are localized 
at the positions of the impurities. Therefore this is just the opposite mechanism 
from that in the Lifshits case. It would be interesting to generalize this approach 
to higher dimensions where similar behaviour can also occur. 

For /U = 0, since the positive and negative S functions play a symmetric role 
this argument doesn't apply anymore. A study of the low energy states on a 
finite interval shows that the existence of quasi zero-modes can account for the 
logarithmic behaviour of the density of states. 

Imposing Dirichlet boundary conditions on a finite interval [—R, R] one finds 
that the ground state energy is (Comtet et al. 1995, Monthus et al. 1996) 

J^rOX Wo\x ) \J-Ri^ Jo / 
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This result is obtained by approximating the true ground state wave function 
near the boundaries by a suitable linear combination of the two linearly inde- 
pendent solutions of -ff+V" = 0: 

^o(a;) = er''^'^(^'' (3) 

and 

dx' 

If (j){x) is a white noise or a random telegraph process with Poissonian lengths, 
the typical behaviour of V'o(a^), given by the central limit theorem, is ipoix) ^ 
e^^. Replacing ipaix) by its typical behaviour one finds that the energy Eq is 
exponentially small in the length of the system Eq ~ e~^. Therefore a quasi 
zero mode of energy E has a typical spatial extension 2R such that R In^ E. 
Coming back to the whole line one finds that the number of such states per unit 
length is 



2R In^E' 

Obviously this argument may be generalized to any one- dimensional disordered 
system for which the zero energy wave functions can be expressed in terms of the 
potential. This is in particular the case of the Anderson model with off-diagonal 
disorder. The discrete Schrodinger equation may be written in the form 

Pn+l^n+l + f3nVn-l = Elfin (5) 

where /3„ are random variables. The model of Dyson of an harmonic chain with 
random masses belongs to this class. For any configuration of the disorder one 
can write down two independent zero energy solutions, One of them is obtained 
by solving the recurrence relation (5) with E = 0. A zero energy state satisfying 
the boundary conditions <po = 1 and = is 

n/32fc-l 

If the P„ are independent identically distributed random variables, the typical 
behaviour of (fin, given by the central limit theorem, is again of the form |</?n| 
Q±Vn_ One will therefore get the same low energy behaviour as before. For earlier 
references see (Theodorou and Cohen 1976, Markos 1988, Bovier 1989); another 
derivation of the logarithmic singularity is given in (Eggarter and Riedinger 
1978). 
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2.3 Remarks 

1. All these arguments are based on typical realizations of disorder. There exist 
however certain quantities whose behaviour cannot be obtained by this type 
of reasoning. This is in particular the case of the average ground state energy. 
Its dependence on the size of the sample has been obtained in (Monthus et 
al. 1996). It is given by the stretched-exponential function (Eq) = exp {—R^) 
where the exponent ^ depends only on the nature of the correlations in the 
potential (the Gaussian white noise corresponds to ^ = 1/3). In (Monthus 
et al. 1996) it is shown that this behaviour is indeed supported by atypical 
realizations of the random potential. 

2. The existence of a singular behaviour in the density of states implies, by the 
Thouless formula, that there will be a corresponding singularity in the local- 
ization length \{E). This quantity indeed diverges as InE which reflects the 
appearance of a critical state at E = 0. It was recently pointed out (Steiner 
et al. 1997, Balents and Fisher 1997) that there exists another length scale in 
the system - the correlation length which controls the decay of the two-point 
Green's function. It behaves like In^ E and thus diverges faster than the 
localization length. These results are in agreement with those presented in 
(Bouchaud et al. 1990); although the full correlations were not computed ex- 
actly, it is shown in this paper that the two-point Green's function is indeed 
given by 

{x\— — -\y) ~ >^c„e 2i„^E I "I 
n — tj E^o '-^ 

n 

It would be interesting to compare this method with Berezinskii's diagram- 
matic technique (Berezinskii 1974) recently used in (Steiner et al. 1997). 

3. The existence of a critical state at S = is best understood when this model 
is reinterpreted in the context of classical diffusion in a random medium. A 
diffusive behaviour at large time requires the existence of an extended state 
at .B = (Tossatti 1990). 

4. The thermodynamic properties of some one-dimensional spin systems with 
random exchange couplings can be reinterpreted by using a mapping of the 
spin system onto a model of free fermions. This approach, which can be 
traced back to the pioneering work of Lieb, Schultz and Mattis (Lieb et al. 
1961), was used by Smith (Smith 1970) in the context of the X-Y model. 
Exact result for quantum phase transition in random X-Y spin chains with 
a comparison with the renormalization group approach (Fisher 1994) were 
obtained by McKenzie (McKenzie 1996). Quite recently, a similar approach 
with a different type of disorder was developed by Fabrizio and Melin (Fab- 
rizio and Melin 1997). It is also worth mentioning the nice paper of Steiner, 
Fabrizio and Gogolin (Steiner et al. 1997) extending this analysis to the case 
of correlations and boundary effects. 
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3 Correlation functions of the ground state wave function 

The localization properties of the wave function can be characterized by the 
density-density correlation function. Various techniques have been developed, 
mainly by the Russian school (Lifshits et al. 1988), to compute such quantities 
in the weak disorder limit. In the supersymmetric model, one may take advantage 
of the fact that the ground state wave function is known exactly as a functional 
of the disorder. If the disordered potential is white noise, this allows one to 
compute the corresponding n-point function by using a mapping with Liouville 
quantum mechanics. Such a calculation was recently carried out by Shelton and 
Tsvelik (Shelton and Tsvelik 1997) for the case n = 2 and n = 3 in the context 
of spin Peierls systems. An extension of this result to arbitrary n is given below. 
We first consider periodic boundary conditions. For complctncss wc also give the 
corresponding formula due to Broderix and Kree in the case of free boundary 
conditions (Broderix and Kree 1995). 



3.1 Periodic boundary conditions 

We consider the supersymmetric Hamiltonian (1) in which wc set ^^^^^-^ = 
Jq dx' 4>{x') for consistency with the notation of previous work (Monthus and 
Comtet 1994, Comtet et al. 1996). We are interested in the following two sections 
in the statistical properties of the zero mode wave function 



jfdx'e0u(^' 



1/2 



(6) 



when the disordered potential 4>{x) is white noise. We consider a system of length 
L and impose periodic boundary conditions. This is achieved if the disordered 
potential is a Brownian bridge (U{Q) = U{L)). The average over the disorder is 
performed through the Wiener measure 



pU{L)=0 J „L /d[/(x)\2 

{■■■)=Af VU{x) ■■■e~^Jo '^''y^) (7) 

Ju(o)=o 

where is a normalization to be determined. Our aim is to compute the corre- 
lation functions of the square wave function iI}q{x): 

Cn{xu---,Xn) = {\Mxi)?---\MXn)?) ■ (8) 

In order to perform the average over U [x) it is convenient to exponentiate the 
denominator coming from the normalization of the wave function. By using the 
integral representation of the F function, one gets 



r{n) 



' dpp" 



VU{x)e [^'^^^ 

U{0)=0 



+0U{xi) + -+pUiXr,) 
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This expression establishes a link with one-dimensional Liouville quantum me- 
chanics. For a detailed discussion of this model we refer to (Monthus and Comtet 
1994) and (Kolokolov 1993, Kolokolov 1994). Considerable simplification occurs 

if one makes use of the fact that, in this theory, a change in the coupling con- 
stant can be interpreted as a translation in U space. This suggests the change 
of variable p = ae^^ , where a = which leads to 

9^/w"n-"+5 /■+°° rU{L)=U 

1 (n) Ju{o)=u 

(9) 

where the action 



-L 



L 

dx 



2a \ dx 



(10) 



is associated with the Liouville Hamiltonian 



The role of the normalization constant M is to insure the relation 

dXi - ■ -dXnCnixi,- ■ ■ ,Xn) = 1- 



Since J\f is independent of n this relation may be used in the case n = 1 to find 
M. 

One may express the correlation function in terms of the Liouville propagator 
G^{U,U') = {U\e-''"^\U'). Choosing L > xi > X2 > ■ ■ ■ > Xn > where 
Xij = Xi — Xj one gets 

C„(a;i,---,ar„) = ^-— / dC/idC/2 • • • d[/„ Gl_,,„(C/„, C/i)e'5^^ 

G^,,{Ui,U2)e'^^' ■ ■ ■ G,„_,„(C/„_i, Un)e^^'' 
Using the eigenstates | ) of the Liouville Hamiltonian^ gives 
1\I'k1 (T^+i /■+°o 

Gn(a:i,---,x„)= . / d/ei---dfc„(/c„|e''^|fci)---(fc„_i|e^^|fc„) 

a(I.-xi„) 2 . 2 

e 4 '^i 4 (12) 

Knowledge of the wave functions gives the matrix elements 

1 , i,2 _ 1.12 

{ k le'^^l k' ) = -Vkk' sinh Trfc sinh Trfc' — - — ■ — — -. (13) 

8 cosh TTK cosh Trfe 



^ The Liouville Hamiltonian has a continuous spectrum Hi,ijik{U) = ^j-ipkiU) where 
the wave function is il>k{U) = V^h^^^j^.^^ ^2c~ 
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Expression (12) allows one to get the long range behaviour (when all the dis- 
tances involved are large compared to i). Using the Laplace method one even- 
tually finds 

^n(»^l5 ' ' ' 1 *^n) — 



n — l ' 

{Ana) — r{n) [{L - xi + .-r„)(.'j;i - X2) ■ ■ ■ {x^-i - a;„)] 

(14) 

Despite the coefficients being different to those found by Shelton and Tsvelik 
(Shelton and Tsvelik 1997) when n = 2 and n = 3, we get the same behaviour 
as a function of the distances. This expression shows the existence of long range 
correlations. A nice interpretation of the algebraic tail was recently given in 
(Laloux and Le Doussal 1997). It is suggested that the exponent 3/2 is associated 
with configurations of the disorder where U (x) returns to its starting point. 
It is interesting to point out that the same exponent appears in conventional 
localization theory. In this case wave functions with the same (Berezinskii 1974, 
Gogolin 1975) or nearly equal (Gor'kov et al. 1983) energy are weakly correlated 
at large distances. The correlation function decays exponentially with the power 

law preexponential factor (-)^^^. 



3.2 Free boundciry conditions 

Computation of correlation functions of the type in (8) was first performed by 

Broderix and Kree (Broderix and Kree 1995) in a different context: |V-'o(a;)P is 
interpreted as the equilibrium Gibbs measure for the potential U{x). Setting 
C/(0) = and letting U{L) be free they get 



C^^{XI,--- ,Xn) ^ — - 

7r(47ra)V- r(n) \/xn{L - xi) [{xi - 0:2) • • • {xn-i - Xn)]^'"^ 

(15) 

which is very similar to (14). The different boundary prescription therefore in- 
duces a slight change in the x dependance near the boundary. 



4 Multifractality 

In a recent series of publications, de Chamon et al. (de Chamon et al. 1996) have 
considered a two-dimensional localization model which exhibits a localization 
transition at zero energy. They consider a two-dimensional Dirac Hamiltonian 
in a random magnetic field -B(r) = A(j){v). In this case the zero energy solution 
(Aharonov and Casher 1979) can be constructed explicitly for any realization of 
4>{t). They consider the particular "groimd state" solution ■^(r) a (e"'^*-''-', 0). For 
a Gaussian disorder P[4>] = cxp — ^ J dr (V(/))^ the successive moments of ip{r) 
are encoded in the partition function 2{q) = / dr e"^''^'''-' since / dr^/;^^(r) = 
■^^^ ■ The multifractal exponents of the wave function can be obtained by using 
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a formal equivalence with the problem of a directed polymer on the Cayley tree 
(Derrida and Spohn 1988). 

In the following we consider a one-dimensional version of this model with, 
however, a different type of disorder and take into account the wave function 
normalization as in (Kogan et al. 1996). The starting point is the two-dimensional 
Euclidean Dirac operator 

11/)= iai {dy +iAy)+ 102 {d^ + iA^) 

(Ji are the Pauli matrices and the gauge field is given by 



{ 



Ay = fix) 
= 0. 



We may take eigenstates of the form 'tjj{x,y) = (x(a;)e"^^, O). The eigenvalue 
equation —^x = ^^X then becomes 

[-dl + + f{x)f - a^fix)] x{x) = E\{x). 

We are thus led to a onc-dimonsional Schrodingcr equation with a supcrsym- 
metric potential V{x) = f'^ix) — a3(p'{x) where ip{x) = uj + f{x). 

In an earlier work by one of us (Comtet ct al. 1988) this approach was used 
to study the density of states of the two-dimensional Dirac operator. If f{x) is 
white noise we may use the density of states of the one-dimensional problem and 
integrate over the free motion on the y axis. The resulting expression displays 
an enhancement at low energy as compared to the free case. 

Here our purpose is to characterize the fluctuation properties of the ground 
state wave function x.{->') = c ^ ^^yf^y\ By rising two different approaches for 
two boundary prescription we compute exactly the successive moments of the 
normalized ground state 

^o(a;) = — ^ 

/o <^x'x^{x') 

In order to keep unified conventions we will parametrize Vo(a;) as in equation 
(6). 



4.1 Moments of xj^Q^x) for periodic boundary conditions 

If the random potential obeys periodic boundary conditions, the situation is the 
one described in section 3. The moment of order 2n is then related to the n-point 
correlation function at coinciding points 



(|Vo(:r)|'")=C„(:r,---,x). 
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In terms of the Liouville propagator this is expressed as 

Tinj 



Jo 



where the matrix element is 



^ ^ m=0 



The moments are eventually given by 



/|^o(x)|2") = -—^1— V a" 2^(2771- l)!!-4^ (16) 
Xiv^uv ;i / 2"(2n-l)!!^ ^ ' (aL)"" ^ ' 

^ ' m=l ^ ' 

where the coefficients aj^ are defined by the equation 

n— 1 n 
m=0 m=l 

For example a" = 1 and a" = -r(n)^. 

One may extract from this expression the asymptotic dependance of the 
moments (when L is large compared to ^) 



4.2 Moments of tpo{x) for free boundciry conditions 

If one leaves U [L) free instead of imposing periodic boundary, then the Brownian 

that enters in the wave function (6) is no longer a Brownian bridge but a free 
Brownian motion starting from U (0) = 0. 

In the average over the disordered potential one must now take into account 
all the Brownian paths starting from without any restriction on the final point. 
Except for this slight modification the formalism is the same as in previous 
sections. One is led to 

(V'o'"(^i)) = ^ dC/ / VU{x) e-^-e"''^(-) 

W J-oo Ju{0)=U 

= ^ / dC/dC/idC/2 Gl-,{U2, Ui)e^^'''G,{Ui,U). 

It is interesting to present another derivation of the moments using the lan- 
guage of previous works (Oshanin et al. 1993a, Monthus and Comtet 1994, 
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Comtet et al. 1996) devoted to the study of the statistical properties of the 
exponential functional Z^^^ = J^^ dx e^(^"^+^^^('^)' . Here \i is the drift for 
the Brownian motion pU{x) = —{ajix + \/2aW{x)) and W{x) is a Brownian 
motion of average and variance 1 defined on [0, L\ . 

We will prove that the moments of '4!q{x) can be expressed in terms of the 
characteristic function 0^^' (p, L) = (c^p^'l ^ ) (Laplace transform of the distri- 
bution law of ^jf')- For this purpose we may note that ^^^{x) involves two 
exponential functionals. We may write 

'4^o{x) = —f 

/o dye/3(^(!')-f^(^)) 

and separate the denominator into two parts 

^yQp{U{y)-U{x)) ^ r ^^yQ-i-tJ.av+V2^B{y)] (Jy g- [May+ V2^B(y)] 

Jo Jo Jo 

where B{y) = W{x — y) — W{x) and B{y) = W{x + y) — W{x), respectively 
defined on y G [0, x] and y G [0, L — x], arc two independent Brownian motions 
starting from zero -B(O) = -B(O) = 0. The denominator is then a sum of two 
statistically independent exponential functionals zi~'^^ and Z^}^. 
The moments may then be rewritten as 

1 \ 




dpp"- V^""^ {P, x) (t>^^^ {p, L-x). 

Where the characteristic functions are given in (Monthus and Comtet 1994, 
Comtet et al. 1996). In the case fi = which is of interest for us 

<^(°)(p,L) = - Tds cosh^ /fi, (2^1^) e-^^\ (18) 
71" Jo 2 V V a/ 

One may extract from this expression the dominant behaviour when ax ^ 1 
and a{L — x)^ 1 

One may check that this result agrees with the one of Broderix and Kree (15) 
in the case n = \. Whenever x belongs to the interval [0,L], but is not on the 
edges, the moments still behave as {tpo"{x)) ~ 7;. 
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We may also explore the behaviour of the moments of the wave function on 
the edges (?/'q"(0)) = (V''o"(^))- This quantity is equal to the moments of the 
partition function zf^ for negative orders 



We may easily extract the dominant behaviour of this quantity with the help of 
(18) and eventually find 

(^o^"(0)) - (20) 

This shows that the wave function fluctuates more on the edges of the interval 
(when x^^oiL — x<^^) than in the bulk (when a; ^ ^ and L — x ^ ^). 

The square of the wave function on the edge is interpreted in the problem of 
classical diffusion in the random potential U{x), as the steady current density 
when a constant density of particle is imposed at a; = in the presence of a trap 
at X = L. The distribution of the steady current was found in (Oshanin et al. 
1993b, Monthus and Comtet 1994). Setting J = ^V'o(O) distribution of J is 
a log-normal law for small values of J 

P(j) ]_ e-^^"'-^ 

and possesses an exponential tail for large values of J 

1 1 ,2 , 
P(J) ^ e^-^ 
J— >oo yiraL J 

which is responsible of the behaviour of the moments given in (20). 
4.3 Discussion of the results 

We have just seen that the one-dimensional model is closely related to the model 
studied in (de Chamon et al. 1996). Since the one-dimensional model is exactly 
solvable, it is interesting to study its multifractal properties, although one cannot 
expect it to give the same behaviour. In fact, in both models the potential is 
long-range correlated; though in (i = 2 the correlations are logarithmic whereas 
in d = 1 they are linear. 

We define the scaling exponent f(g) which characterizes the behaviour of the 
critical wave function at = as 

L — ^OC' 

Using equation (17) we get f(g) = 0. This exponent is the scaling exponent of 
the average moments of the wave function. We may introduce a slightly different 
scaling exponent 

L^oo ln(l/L) 
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which is the mean exponent of the critical wave function. In terms of the partition 
function z'f^ {(3) = dx e^^^^^ introduced before, it may be rewritten 



r(g) = lim 



InZf (9/3))-g(lnZf 



L^oo ln(l/Z/) 

The mean free energy that appears has been calculated in (Comtet et al. 1996) 
to be 




lnZ^°\f3)) = 2J— + C -\na ^ + 



1 



where C = — r"(l) is the Euler-Mascheroni constant. Again, one eventually 
gets a scaling exponent T{q) = which agrees with the behaviour found in (de 
Chamon et al. 1996) in the strong disorder regime. By using the mapping with 
the random directed polymer model (Derrida and Spohn 1988), de Chamon et 
al. have shown that this regime corresponds in fact to the low temperature phase 
of this model. The existence of such a link also appears in the one-dimensional 
case. A comparison of the moments of Z^^ {j3) with those of the random energy 
model (REM) (Derrida 1981) reveals some striking similarities. For the one- 
dimensional case the expression of the moments given in (Monthus and Comtet 
1994) read 

/i'70\n\_ 1 i -^("-) Y^/ ixn-fe^aLfe^^fe+n . 1 

In particular, for large L one obtains 

These moments grow in the same manner as in the REM. The main difference is 
that, unlike in the REM^ there is here no transition above which the behaviour 
would change from e"^" to e"'^". In some sense one can consider that the tran- 
sition temperature is sent to infinity. This explains why, in this one-dimensional 
case, one only probes the low temperature phase. It would be extremely interest- 
ing to explore intermediate cases with weaker correlations, as recently suggested 
in (Bouchaud and Mezard 1997). 
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